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A SPIN MODULATED TELESCOPE TO MAKE TWO DIMENSIONAL CMB MAPS 



ON 

o 

Q 



> 

(N 

(N 
(N 

On 
On 



Oh' 
i 1 

O 

CO 



X 
3 



J. Staren, P. Meinhold, J. Childers, M. Lim 1 , A. Levy, P. Lubin 

Physics Department, University of California at Santa Barbara, Santa Barbara, CA 93106-2530 

M. Seiffert, T. Gaier 

Jet Propulsion Laboratories 
N. FlGUEIREDO 

Escola Federal de Engenharia de Itajuba, 37500-000, Itajuba, MG, Brazil 
T. VlLLELA, C.A. WUENSCHE 
Divisao de Astrofi'sica, Institute National de Pesquisas Espaciais, 12225-010, Sao Jose dos Campos, SP, Brazil 

M. Tegmark and A. de Oliveira-Costa 

Department of Physics, University of Pennsylvania, Philadelphia, PA 19104 
Draft version February 1, 2008 

ABSTRACT 

We describe the HEMT Advanced Cosmic Microwave Explorer (HACME) , a balloon borne experiment 
designed to measure sub-degree scale Cosmic Microwave Background anisotropy over hundreds of square 
degrees, using a unique two dimensional scanning strategy. A spinning flat mirror that is canted relative 
to its spin axis modulates the direction of beam response in a nearly elliptical path on the sky. The 
experiment was successfully flown in February of 1996, achieving near laboratory performance for several 
hours at float altitude. A map free of instrumental systematic effects is produced for a 3.5 hour observa- 
tion of 630 square degrees, resulting in a flat band power upper limit of (£(£+ l)Cg/ (2ir)) - 5 < 77 /iK at 
I = 38^20 (95% confidence). The experiment design, flight operations and data, including atmospheric 
effects and noise performance, are discussed. 

Subject headings: cosmic microwave background - cosmology: observations - instrumentation: 
detectors 



1. INTRODUCTION 



The Cosmic Microwave Background (CMB) is a unique 
relic of the primordial Universe and the study of its prop- 
erties promises to advance our understanding of the for- 
mation of structure in the universe. Detection of CMB 
anisotropy helps to discriminate among potential models 



tion over a range of microwave f r equencies (Coble et al 
1999|,|de Oliveira-Costa et al. 1998|, [Netterfield et al. 1997|) 



with different cosmological parameters (Bond, Efstathiou 



fc Tegmark 1997 ; Imigman et al. 1996 ; Hu ; Sugiyamg 



fc Silk 1997| ; [Zaldarriaga, Spcrgcl fc~SeTjak 1997| ). Much 
of the discriminating power of CMB measurements comes 
from angular scales below 1 degree, or Legendre multi- 
poles, I, greater than 100 in the angular power spectrum. 

The effort to understand the power spectrum of the 
CMB for large I has c ontinued since the discovery of 



10021, I Schuetor et al. 10051, I Roadhoad ot al. lOSO E 



aniao tropy in the CMB ( gmoot et al. — 1992 ;- A kiop et 



-al 



Gmidor 



sen ct al. 19951, |Lim et al 1996| , |Cheng et al. 1997 



Netter 



With sensitivity from the beamsize to the large angle of 
the spinning chop, a multi-pixel spin-modulated telescope 
easily covers enough sky area to make precision measure- 
ments of the angular power spectrum. In principle, future 
experiments will enable cosmological parameter extraction 
from such power spectrum measurements. 

The HACME experiment is the precursor to several 
spin modulated experiments that will accurately measure 
the angular power spectrum using cryogenic HEMT (High 
Electon-Mobility Transistor) amplifier receivers flown on 
balloon platforms to minimize atmospheric contamination 
and enable large angle chopping. HACME has a single 
Q-band r eceiver fed by the ACME telescope (Meinhold et 



field ct al. 1997). Many previous experiments minimized 



errors for the best sensitivity at a single angular scale on 
a small region of sky. An ideal experiment will have suf- 
ficient sky coverage with enough sensitivity over a broad 
angular range to measure the power spectrum at several 
angular scales. Broadly speaking, the optimal strategy 
is to aim for a signal to noise of order unity per pixel 
with the largest and most symmetric sky coverage possi- 



blc, f or a given resolution and instrument sensitivity (Knox 



1997 



). Instrumentation has evolved so that sensitive 
large area maps can be obtained with sub-degree resolu- 



al. 1992) and modulated by a 1.2 m rotating canted flat 
mirror. The single horn and conventional duration bal- 
loon flight allows an upper limit on CMB anisotropy and 
upcoming instruments will have more horns and longer 
flights to make power spectrum estimates. The novel op- 
tical spinning technique, described in detail in §2.1, allows 
rapid two dimensional modulation of the instrument to 
reduce the effects of 1 /f noise in the instrument while ty- 
ing together widely separated regions of the sky in a short 
time. Changes in optics, computers and data acquisition 
made for HACME are also described in §2. In §3, the 
observations made and the data collected are discussed. 

2. INSTRUMENT 
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2.1. Optics and Scan Strategy 

HACME's optical system (Figure [[]) consists of a fixed 
ellipsoidal mirror, the original ACME 1 m parabolic pri- 
mary mirror, a new 1.2 m spinning flat mirror and a baf- 
fle surrounding the flat mirror between 2° to 5° from the 
optical axes to intercept direct spillover paths from the 
Earth to the edge of the primary mirror. The spinning 
flat mirror, canted 2.5° relative to its spin axis, modulates 
the telescope beam in a nearly elliptical path on the sky 
so we become sensitive to anisotropy on scales from the 
beamsize, 0.8° FWHM, to the size of the major axis of the 
ellipse, 10.0°. The polar angle, 8 r , and an azimuthal angle 
about the center of the path, <f> r , are 



s(0 r ) = cos(26» c ) + 2sin 2 (6»/)! 



\<f>f) (1) 



and 



tan( 



sinO/) 



cos(6> c ) cos(0j) — sin(6* c ) sin(0j) cos(</>/) 



cos(# c ) cos(</>/) + | sin(20r) sin(# c ) sm 2 (4>f) 

' (2) 

where 8 C = 2.5° is the cant, 9i — 31.9° is the angle be- 
tween the mirror spin axis and the main optical axis and 
(fif is the phase of the flat mirror spin. Performing an 
azimuth scan adds another level of modulation, building 
up a two dimensional image of the CMB. At the start of 
the azimuth scan, the leftmost part of the beam path is 
aligned with a fiducial star. The gondola scans 20° to the 
right and returns at a rate of 0.17 degrees per second. This 
process continues as we track the fiducial star across the 
sky. This observing strategy is shown on a map in Figure 
|^. Data taken in this manner produce l arge maps w ith 
nearly diagonal noise covariance matrices ( Wright 1996| , ) . 
The resulting region of sky covered is almost 30° wide and 
between 10° and 30° high for a typical three to five hour 
observation. 




FIG. [l]. Schematic Drawing of the Major Optical Elements 
of HACME: A. Baffle; B. Canted Spinning Flat Mirror; C. 1 m 
Parabolic Mirror; D. Ellipsoidal Mirror; E. Dewar; F. CCD Cam- 



Telescope response is measured on an antenna range in 
Santa Barbara with a Gunn oscillator source. The best fit 
Gaussian to the main lobe along the azimuthal and ele- 
vation axes at 40 GHz are 0.84 ± 0.02° and 0.75 ± 0.02° 



FWHM respectively. The beam shape is independent of 
the orientation of the flat mirror to within the measure- 
ment errors and varies as the inverse of the frequency of 
the subbands described below. Observations of the Moon 
during flight confirm our ground tests of the beam size. 
We measured the illumination taper of the primary to be 
below -25 dB around the edge and calculated the taper 
around the rotating flat mirror to be always less than - 
39 dB. Sidelobe measurements made with the main baf- 
fle in place demonstrate response below -60 dB for angles 
greater than 20° from boresight along the elevation axis. 

2.2. Receiver and Data Acquisition 

The HACME Q-band rec eiver is the same as that of 
SP94 ( pundcrscn ct al. 1995 ), with a five-stage InP-GaAs 
HEMT multiplexed into three 2.3 GHz bandwidth sub- 
bands centered at 39.15, 41.45 and 43.75 GHz, but is used 
as a total-power radiometer for this instrument. The out- 
put diodes for each subband connect to a low noise pre- 
amplifier, an A/C coupling stage (1/RC w 7 sec) and a 
linear voltage to frequency (V/F) converter. The output 
of the V/F is counted synchronously with the spin of the 
flat mirror, 125 triggers per rotation. An indexed opti- 
cal encoder attached to the spinning mirror shaft provides 
the triggers that divide each mirror spin into 125 "sectors" . 
The mirror spins at 2.500 ± 0.004 Hz, resulting in an ideal 
integration of the diode outputs at 3.2 msec. Pulses from 
a 16 MHz clock are counted for normalization. Pointing 
and receiver data are stored on board the payload and a 
compressed form of the data is sent to the ground. Analy- 
sis is performed on the data recovered from the hard disk 
after landing. 

2.3. Calibration 

The telescope is calibrated using two methods. First, 
the absolute calibration, including atmospheric attenua- 
tion and telescope loss, is determined by observing the 
Moon. A model of lunar emission by S. Keihm (JPL) is 
compared to two-dimensional maps made from the data. 
Secondly, an ambient temperature "Eccosorb" target can 
intercept the beam between the primary and secondary 
mirrors to provide relative calibration during the flight. 
The latter calibrations are performed about every hour 
during flight and determine the relative calibration to 
±3%. The following results are calibrated in antenna tem- 
perature units. Conversion of the data to CMB thermo- 
dynamic temperature requires multiplication by 1.04, 1.05 
and 1.05 for the 39, 41 and 43 GHz subbands. 

3. DATA 

3.1. Observations 

This experiment flew for the first time on 1996 Febru- 
ary 11 from Ft Sumner NM, obtaining 15 hours at a float 
altitude of 35.4 km. Operations proceeded smoothly, with 
detectors and all servo systems reaching expected perfor- 
mance for the entire flight. Of the 12 nighttime hours at 
float, 8 hours were spent mapping three extended regions 
near the stars Alpha Leonis, Gamma Eridani and Gamma 
Ursae Minoris. These regions were chosen for low fore- 
ground emission and with systematic tests in mind. All 
regions observed are between 40° and 50° from the galac- 
tic plane. When observed, aLeo was rising in the east, 
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7Eri was setting in the southwest and 7UMi was transit- 
ing in the north. The remainder of the time was spent 
on calibrations, main lobe mapping using the Moon, sys- 
tematic tests such as operating without the servo systems 
and RF transmitters. A second flight was attempted on 
1996 June 1 but was terminated by a balloon failure before 
reaching float altitude. 
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FIG. Gamma Ursae Minoris map. The map is made from the 
43 GHz subband data. There are over 3000 gnomonic 28' square 
pixels. The contours encompass regions covered greater than 12.8, 
6.4 and 3.2 seconds of integration per square degree. The white dots 
represent the 125 samples taken along an elliptical path for each spin 
of the flat mirror. Spins at each end of the azimuth scan are shown 
along with the dashed boundary of the region observed during one 
scan of the telescope. 



spin-synchronous offset for each subband. The peak-to- 
peak amplitudes of the offsets are 1.3, 1.6 and 4.0 mK in 
the 39, 41 and 43 GHz subbands, respectively. The ampli- 
tude of these signals varies < 5%/hr over the entire flight. 
Among the possible causes are atmospheric emission, ther- 
mal emission from the optics, the CMB dipole and far 
sidelobe contamination. Because the observed offset varies 
slowly, we are able to remove it on short timescales. We 
discuss the possible contributions to this offset and its re- 
moval below. 
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FIG. ^. Spin-sychronous offset in 43 GHz subband vs. Altitude. 
The diamonds are the zenith temperatures assuming the detected 
offset is due only to the atmosphere. The dotted line is an atmo- 
spheric model prediction using Zhevakin-Naumov line shapes. The 
dashed lines is for Van Vleck-Weisskopf lines shapes. The observa- 
tion and models differ at float altitude because the atmosphere is no 
longer the dominant contribution to the offset. 



3.2. Flight Performance 

During flight, the HACME telescope performed as ex- 
pected with no known difficulties. Once scanning a region, 
the telescope collects data very efficiently. For the 7UMI 
observation, after removal of 16% of the data due to cali- 
brations, we are able to use 99% of the remainder. One of 
the advantages of a balloon environment is that essentially 
all of the data is useable as there is no "bad weather" . The 
power spectra of noise for the three subbands have high 
frequency (white noise) limits of 0.9, 0.5 and 0.4 mK v / sec 
with a typical knee frequency of 60 Hz (the frequency at 
which the noise power is double the white noise limit). 
The subbands have typical correlations between 0.5 and 
0.7 due to HEMT amplifier gain fluctuations. This perfor- 
mance is consistent with our laboratory measurements of 
this receiver and is typical of cryogenic HEMT amplifiers. 

3.3. Systematics 

3.3.1. Synchronous Offset 

In order to make maps from the data, we need to re- 
move systematic offsets. When the data are binned into 
"sectors" rather than sky pixels, we find an anomolous 



While the atmosphere is the dominant cause of a spin- 
synchronous offset on the ground, we find that it is small 
at our float altitude and remarkably consistent with our 
predictions at lower altitudes based on a model for atmo- 
spheric emission. We verify the spectrum of the atmo- 
spheric contribution and its elevation and altitude depen- 
dence up to 25 km with data taken during the ascent of the 
balloon payload (Figure ||). At an altitude of 35 km, the 
largest peak-to-peak signals predicted are 500, 600 and 800 
/iK for our detection subbands. During the aLeo measure- 
ment, the telescope tracked a 22° range in elevation, over 
which the atmospheric signal would have doubled. With 
no detectable monotonic change in the offset, the atmo- 
spheric signal in the 43 GHz band is limited to less than 
600 /iK for this observation. 

Thermal emission from the optics is expected to be a 
larger contribution to the spin-synchronous offset and op- 
posite in phase from the atmosphere. Our Q-band receiver 
detects only one linear polarization and a synchronous sig- 
nal is expected as the th e angle of the flat mir ror relative 
to the receiver changes (Wollack et al. 1997). For our 
horizontal E-plane, the signal is largest when the beam 
elevation is greatest. This signal depends upon the tem- 
perature and emissivity of the flat mirror. The peak-to- 
peak amplitude of this signal is predicted to be 7.4, 7.6 
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and 7.8 mK in the increasing frequency subbands. The 
polarization signal is calculated to decrease just 10% over 
the course of the entire flight as the flat mirror cools. 
We calculate the signal expected from the motion- 



induced doppler dipole anisotropy (Lineweaver et al 
19960 For a 10° throw, the maximum possible signal from 
the dipole measured by HACME would be 580/jK. The 
synchronous offset produced by the dipole is equal in the 
three subbands and its phase and magnitude varies with 
the equatorial coodinates of the target. 

A linear combination of the sources above cannot ex- 
plain the observed offsets throughout the flight however, 
so other effects must be present. Far sidelobe contamina- 
tion, Radio Frequency Interference, chopped spillover and 
electrical interference are possibilities. Because the Earth 
makes the largest contribution to far sidelobe response, we 
expect this systematic error to be elevation dependent. Be- 
cause the spin-synchronous offset did not change apprecia- 
bly during the aLeo observation, sidelobe contamination 
is not likely a major cause of systematic error. Radio Fre- 
quency Interference was of great concern to us, due to the 
necessity of the NSBF flight transmitters and our previous 
experience. We were however able to turn off and on all 
on board transmitters (we even obtained several minutes 
with all transmitters off), and saw no effect on the data. 
We cannot entirely rule out chopped spillover and electri- 
cal pickup. Chopped spillover refers to near field off-axis 
telescope response modulated by the flat mirror or its bal- 
ance arm. We performed careful noise integration tests 
prior to flight to levels far below the observed offset, but 
differences between ground tests and flight could lead to 
systematic errors from electrical interference. 

3.3.2. Offset Removal 

Current data does not allow the separation and iden- 
tification of the modeled contributions partly due to the 
cancellation of the effects. A conservative approach to off- 
set removal is taken because the offsets detected are due to 
a superposition of the contributions listed above. We re- 
move all data at the first and second harmonic of the spin 
rate before analysis, as the offset is expected and observed 
to be significant only at these harmonics. Only 3.2% of 
the 2 to 3 million degrees of freedom in an observation are 
sacrificed to this harmonic removal. We find it a robust 
approach that sacrifices sensitivity to multipoles £ < 20 
and use it for all following analysis. 

3.4. Maps 

The total power data collected by HACME are easily 
converted into sky maps. Measurements of each "sec- 
tor" are assigned a sky position using Eqns. [l] and || and 



the azimuth and elevation of the center of each flat mir- 
ror spin. Averaging measurements of sky pixels results 
in maps dominated by stripes due to 1// noise and by 
the synchron ous offset. The map- making process in and 
Wright 1996] is needed to reduce striping and excess noise 
due to the amplifier gain fluctuations and is convenient for 
the offset removal. This algorithm accounts for the cor- 
relations introduced by 1// noise and optimally reduces 
striping in the maps by using the inter-connections of our 
scan pattern. The same pointing information is used along 
with 1// noise estimates to make improved maps. Wiener 
filtering using signal-to-noise eigenmodes can be used to 
project out the map modes with the best sensitivity. 
Application of this method to these data is described in 



Tegmark et al. (1999, hereafter T99). Figure H shows a 



Wiener-filtered map from T99 of the 630 square degree, 
3.5 hr, 7UMi observation, with integration time contours 
and the scan pattern overlayed. All structure seen in the 
map is consistent with instrument noise with pixel errors 
a factor of 2.6 greater than those expected for the white 
noise levels indicated above. The sensitivity is best in the 
inner, more heavily sampled regions. T99 find a flat band 
power upper limit of + 1)GV(2tt))^ 5 < 77 (95% 
confidence) at I = 38^20 from this observation. 

4. CONCLUSION 

We have built, tested and flown a novel experiment to 
make large sky area two dimensional maps of the CMB. 
The new equipment used at balloon altitudes, the HEMT 
amplifier and canted spinning flat mirror, worked flaw- 
lessly. We find that atmospheric emission should not sig- 
nificantly contaminate data of future experiments. The 
maps produced yield an upper limit to CMB anisotropy 
because of the poorer sensitivity than expected. They 
show no systematic effects or striping from the 1/f noise 
and demonstrate that the scan stategy used will be an 
excellent way to make more CMB anisotropy maps. 
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